Electric Dipole Spin Resonance for Heavy Holes in Quantum Dots 
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We propose and analyze a new method for manipulation of a heavy hole spin in a quantum 
dot. Due to spin-orbit coupling between states with different orbital momenta and opposite spin 
orientations, an applied rf electric field induces transitions between spin-up and spin-down states. 
This scheme can be used for detection of heavy-hole spin resonance signals, for the control of the 
spin dynamics in two-dimensional systems, and for determining important parameters of heavy-holes 
such as the effective g-factor, mass, spin-orbit coupling constants, spin relaxation and decoherence 
times. 
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Spintronics, or spin-based electronics, is one of the fas- 
cinating and rapidly growing areas in solid state physics 
and modern technology 1]. Exploiting both the charge 
and spin degrees of freedom of carriers, it offers wide 
opportunities for developing devices with unique func- 
tionalities. Operation of single charge or spin is the ul- 
timate limit for such devices. Quantum dots (QD) have 
proven useful for this goal [2j , since experiments on read- 
out and coherent manipulation of a single spin in QDs 
have already been performed successfully Fur- 
thermore, due to supression of the spin-orbit interaction 
(SOI) in QDs 0, la El the electron spin in a QD has 
long relaxation times T\ (up to hundreds of milliseconds) 
lj- IHJ ■ The spin thus is an attractive candidate as 
carrier of quantum information [2|, if its state stays co- 
herent over sufficiently long times (described by the spin 
decoherence time T2). At low magnetic fields, in III- 
V semiconductor QDs there is a rather strong hyperfine 
interaction between an electron spin and surrounding nu- 
clear spins lea ding to a significant degradation of the spin 
coherence |ljl [l2l |lj| . Very recently, based on the idea of 
state narrowing via projective measurements several 
proposals have been made to limit the hyperfine-induced 
decoherence ^3^3; an d, therefore, to increase T2 times 
for electron spins in QDs, which currently range up to /xs 



Electron spin resonance (ESR) provides a powerful tool 
for coherent manipulation of spins, which is of great im- 
portance for spintronics [lfj • By applying short resonant 
microwave pulses, an arbitrary superposition of spin-up 
and spin-down states is created. Rabi oscillations and 
spin-echo experiments are based on this approach |5| . In 
such experiments, the ESR signal can be detected di- 
rectly by measuring the absorption of radio- frequency 
(rf ) power via charge transport through a QD [fl, 
or by u sing optical detection of magnetic resonance tech- 
niques [13, UJi- Usually, ESR methods involve mag- 
netic dipole transitions induced by an oscillating mag- 
netic field. However, there has been a strong revival of 
interest in electric dipole spin resonance (EDSR) con- 



trolled by alternating electric fields 0,111 m 00113 

that provides the ability to manipulate and detect elec- 
tron spins at the nanometer scale and may be useful for 
distinguishing different SOI mechanisms. 

Recently, the idea to use heavy-hole (HH) spins (in- 
stead of electron ones) as carriers of quantum informa- 
tion has generated a lot of interest. On one hand, the 
hyperfine interaction of holes with lattice nuclei is sup- 
pressed, since the valence band has p symmetry. On the 
another hand, the spin relaxation and decoherence for 
holes due to HH-phonon interaction could be compara- 
ble to or even longer than that for electrons in flat QDs 

In this Letter, we propose a new method to study and 
control the spin dynamics of HHs in QDs via EDSR. We 
show that, due to SOI between states with different or- 
bital momenta and opposite spin orientations, an applied 
rf electric field may induce transitions between spin-up 
and spin-down states and can be used for detection of a 
heavy- hole spin resonance signal and manipulation of HH 
spins in 2D systems. In contrast, rf magnetic fields do 
not affect the HH spin dynamics (in dipole approxima- 
tion) and, therefore, are inefficient for manipulating and 
observing spin resonance and Rabi oscillations in HH 2D 
systems. 

We consider a single HH confined in a III-V semicon- 
ductor QD exposed to a static magnetic field B of arbi- 
trary direction with respect to the [001] growth direction 
of a QD. Due to confinement along the growth direction, 
HH and LH subbands are split and, for thin enough QDs, 
SOI between LH and HH states can be taken into ac- 
count perturbatively and, therefore, HHs can be treated 
as quasiparticles with pseudospin S aligned along the 
growth direction (S x — S y — and S z = 3er 2 /2, where 
a z is the Pauli matrix) and with the following effective 
Hamiltonian |26j |. 
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H = —(P^ + P') + U(x,y) + H so --g ±f M B B ± a z . (1) 
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m is the effective 



Here 
|e|A(r)/c, A(r) 



HH mass, P = 
(- y B ± /2,xB ± /2,yB x - xB y 
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(B x ,B y ,B±), U(x,y) is the lateral confinement poten- 
tial of a QD, g± is the component of the g-factor tensor 
along the growth direction, and 

H so = l3P-P+P-(J + + iaPlcj + + 37o ^ B B-Pl(j+ + h.c. 

mo A 

(2) 

is the SOI term of HHs consisting of three contributions: 
the Dresselhaus term (/3) j2|| , the Rashba term (a) (2^ , 
and the last term (due to SOI between LHs and HHs) 
combines two effects: orbital coupling via non-diagonal 
elements in the Luttinger-Kohn Hamiltonian (oc P±) 
and magnetic coupling via non-diagonal elements in the 
Zeeman term (oc B±) 28]. Here, a± — (a x ± i<7 y )/2, 
P± = P x ± iP y , B+_= B x ± iB y , 70 and n are the Lut- 
tinger parameters [28j . too is the free electron mass, A 
is the splitting between LH and HH subbands, and the 
SOI constants (a and P) depend on band parameters, 
confinement along the growth direction, and the energy 
spacing between LH and HH subbands A |26| . Note that 
we neglect the Zeeman splitting due to an in-plane mag- 
netic field B || = (B x ,B y ) (b ecause of strong anisotropy 
of HH g-factor: g<< <§; g± : 2 9j] ) as well as the orbital effect 
of the in-plane magnetic field, which is negligibly small 
for B\\ <C ch/eh 2 |30|, where h is the QD height. 

The Schrodinger equation for Hamiltonian Q without 
SOI term (/fso = 0) and with parabolic lateral confine- 
ment [U(x,y) = moj^x 2 +y 2 )/2] has an exact solution 

mi: 

En^m.s = hfl(n + 1) + hL0 c m/2 — Hujzs/3, (3) 
\n,m,s) = c„, m e— (r/0 |m| e- 2 / 2 ' 2 4:i |m|)/2 (r 2 A 2 )|^) 

where (r, ip) are the polar coordinates, n = 0,1,2,..., 
the magnetic quantum number m = — n, — n — 
2, . . . , 0(1), . . . , n — 2,n (which has the same parity as 
n), s = ±3/2, n = y/uj 2 + uj 2 /4, u j z = g±n B B±/h 
is the Zeeman frequency, I — y^h/mil, c n . m = 
y/{{n - \m\) /2)!/?r ((n + \m\) /2)\/l, and L™(x) is the 
generalized Laguerre polynomial. 

Taking the SOI into account as a perturbation, we find 
that in first-order perturbation theory the two states cor- 
responding to the Zeeman-split ground state level can be 
written as follows: 

|±> = |0,0,±3/2)+i/^=|l,±l,T3/2) + /^=|3,±l,T3/2) 
+a± |3, ±3, T3/2) + 7±B ± |2, ±2, T 3/2) , (5) 

where @f = /3(ml) 3 u±(uj 2 _ + Lo^/fiw^, 7 ± = 
S^/^joK^Bimrfuj-frnQAhu)^ , and ui± = Cl±w c /2 (w D , 
u)^, a ± , and (3^ are not relevant for the present dis- 
cussion). From Eq. J5J, it can be seen that i?so leads to 
coupling of the two lowest states |0, 0, ±3/2) to the states 
with the opposite spin orientations and different orbital 
momenta. Due to this spin-orbit mixing of the HH states, 
the transitions between the states |±) with emission or 
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FIG. 1: (color online). Heavy hole spin relaxation rate 1/Ti 
in a GaAs QD versus an applied perpendicular magnetic field 
B± (the height of a QD is h = 5 nm, the lateral size lo = 
s/h/mujo = 40 nm, k = 1.2, 70 = 2.5, 01 = 2.5 [2<|, and 
the other parameters are given in Ref. 26]). Inset: Energy 
differences of lowest excited levels with respect to the ground 
state -E ,o,+3/2- 

absorption of an acoustic phonon become possible and 
this is the main source of the spin relaxation and de- 
coherence of HHs (because the hyperfine interaction of 
holes with nuclei is suppressed) [26| . The spin relaxation 
between the states |±) due to Rashba spin-orbit (RSO) 
and Dresselhaus spin-orbit (DSO) terms is considered in 
detail in Ref. |26|. Similarly, we can find the contribu- 
tion l/r| of the third SOI term in Eq. J2J to the spin 
relaxation rate 1/Ti: 

(6) 

(see Ref. |26( for details). Note that in contrast to elec- 
trons there are no interference effects between DSO, 
RSO, and SOI due to LH-HH coupling and in-plane mag- 
netic field, therefore, the total spin relaxation rate 1/Ti 
is the sum of rates 1/Ti = l/T^ so + l/T^ so + 1/t|. 
Note that in the limit of low temperatures (ksT <C huj-z,), 
T 2 = 2Ti 0. 

In Fig.^the total spin relaxation rate 1/Ti is plotted as 
a function of perpendicular magnetic field B± . As can be 
seen from this figure and Eq. ©, the SOI due to in-plane 
magnetic field leads to an additional peak which is due to 
anticrossing (2(||33| between the energy levels E +z / 2 
and -E 2j _2,-3/2 ( see the second avoided crossing in the 
Inset). In the Inset, the first (third) avoided crossing 
resulting from DSO (RSO) coupling corresponds to the 
first (third) peak of the spin relaxation curve in Fig. ^ 

Now let us consider methods for the manipulation 
and detection of the HH spin in QDs. For electrons in 
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2D structures, an applied oscillating in-plane magnetic 
field couples spin-up and spin-down states via magnetic- 
dipole transitions and it is commonly used in ESR, 
Rabi oscillation, and spin echo experiments From 
Eq. JSJ it is easy to tell that magnetic-dipole transitions 
(An = 0, Am = 0, and As = ±1) arc forbidden and, 
due to spin-orbit mixing of the states |0, 0, ±3/2) with 
z/?f"|l,±l,=F3/2), electric-dipole transitions (An = ±1, 
Am = ±1, and As = 0) are most likely to occur. There- 
fore, the HHs arc affected by the oscillating electric field 
component and not by the magnetic one. 

We consider a circularly polarized electric field ro- 
tating in the XY-plane with frequency u>: E(t) = 
E(sin uit, — cosut, 0). Therefore, the interaction of HHs 
with the electric field is described by the Hamiltonian 
H E (t) — (\e\E/mui)(cosujtP x + sinutPy). The coupling 
between the states |±) is given by (+\H E (t)\-) = Hf_ = 
(#_+)* = dsoEe- iu,t , where 



dso = (|e|Z/2w)(/3fa;_+#w + ) 



(7) 



is an effective dipole moment of a HH depending on DSO 
coupling constants, perpendicular magnetic field B±, lat- 
eral size of a QD, and frequency of an rf electric field (fif 
and ui± are defined below Eq. JSJ). In the framework of 
the Bloch-Redfield theory |3Jj, the master equation for 
the density matrix in the rotating frame can be written 
as follows: 



p z = 2(d S oE/h)p- - (p z 
P+ = <W- - P+/T 2 , 
P- = S I ip+-2(dsoE/h)p z - p-/T 2 , 



(8) 



. Pi = (W+- 



where p z = p ++ — p--, p+ = p~\ e + 

^( i0+ _ e M^ - p-+e-^), p nm = (n\p\m) 

W |_)Ti is the equilibrium value of p z without rf field, 

W nm is the transition rate from state m to state n, T\ = 
1/(W+- + W-+) is the spin relaxation time, T 2 = 1T\ is 
the spin decoherence time |2^ , d r f = ujz — 5lo — u> is the 
detuning of the rf field, and Sui is the resonance shift due 
to decoherence 0- 

The coupling energy between a HH and an oscillating 
electric field is given by 



(H E (t))=Tr(pH E (t)) 



iso 



E(i), 



(9) 



where d S o = d so (ip-+ - ip+-,p+- + P-+,0) is the 
dipole moment of a HH. Therefore, the rf power P = 
—d(H E (t))/dt = —udsoEp- absorbe d by a HH spin sys- 
tem in a stationary state is given by |3£ 



2uj{d so E) 2 T 2 p T z /h 
1 + <5 r 2 f T 2 2 + {2d so E/h) 2 T 1 T 2 ' 



(10) 



In Fig.0 the dependence of P on a perpendicular mag- 
netic field B± and frequency u> of the oscillating electric 
field is plotted. The rf power P absorbed by the system 
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FIG. 2: (color online). Absorbed power P(meV/s) as a func- 
tion of perpendicular magnetic field B± and rf frequency u> 
(T 2 = 2Ti, E = 2.5V/cm, By = IT, and the other parameters 
are the same as those in Fig. . 



has three resonances and one resonant dip. The first res- 
onance occurs at zero detuning <5 r f = 0. For the QD con- 
sidered, the resonance shift 8lo is negligible {5uj » 1kHz), 
therefore, the first resonance appears when the energy of 
rf radiation equals the Zccman energy of HHs: B 1 ^ 1 = 
fvjj/g±PB- The shape of this resonance (at certain u>) is 
given by P w hwp T j2h\\ + h 2 S^T 2 /{2d so E) 2 T 1 }. At low 
magnetic fields (J5 < Bj 2 ), Slu is the Lamb shift induced 
by zero-point fluctuations of acoustic phonons pjrj l3?l ] 



5u> 



(11) 



The field dependence of 5u> shows up in an effective g±- 
factor which scales approximately as B 1 / 2 . 

If the first and second resonances are well separated 
(lu <C lo- ) , then the absorbed power can be estimated as 



P^2Lo(d so E) 2 p 1 z /(h5 2 { T 2 ) 1 



(12) 



in the region of the second and third resonances and 
the resonant dip. The second resonance corresponds 
to anticrossing of the levels E 0< 0,-3/2 an d -Ei. -1,3/2 ( see 
the first avoided crossing in Inset of Fig. ^| at w_ = 
w z H(| (at B 1 ^ 2 = Huj /g±pB\/l + 2m /g±m). At 
the anticrossing point, there is strong mixing of the 
spin-up and spin-down states and the dipole moment 
of a HH spin system is maximal dgQ X = \e\lu>z/2tu. 
Therefore, the height of the second resonance is given 
by (elivzE) 2 /2hwS 2 { T 2 . The resonant dip appears at 
Bj_ = (hujo/2g±pB)\/2mo/g±m, which corresponds to 
0i W- + Pt^+ = and to zero dipole moment (see 
Eq. The third resonance reflects the peak in the spin 
decoherence rate T 2 _1 due to applied in-plane magnetic 
field (see Fig. at the second anticrossing point (the 
second avoided crossing in Inset of Fig. ^) at 2cj_ = ujz 
(B^ = AhuJo/g±pB\/^ + 4mo/<?±m). From the posi- 
tions of the resonances we can determine g±, m, and 
wo, from the shape and the height of those we can ex- 
tract information about the SOI constants a, (3, and SOI 
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FIG. 3: (color online) . Rabi oscillations at three different val- 
ues of the perpendicular magnetic field: B± = 0.8 T (damped 
fast oscillations), B ± = 0.865 T (dotted line), and B± = 0.5 T 
(solid line). By = 0, 5 r f — 0, E = 1.5 V/cm, and the other 
parameters are the same as those in Fig. 



strength due to in-plane magnetic field (which is pro- 
portional to 7o«;/A). Moreover, we can determine the 
dependence of the spin relaxation and decoherence times 
on B±. 

We now consider dynamics of the HH spin (S z ) = 
(3 /2)Tr (per z ) = (3/2)p z . The Rabi oscillations at zero 
detuning (<5 r f = 0) of the spin are given by 

( Sz ) = 5j + e-^ 1 +T 2 - 1 W2|g_ 5 ^ cosa;Ri 



(d so EfT 2 rp T-f 1 - T- 1 



sinw R i H3) 



+ 



where w R = \J{d so E/h) 2 - (Tf 1 - T 2 -1 ) 2 /4 is the Rabi 
frequency, and = (3/2)p^/[l + {d so E /hfT^]. 

As can be seen from Fig. the spin dynamics is es- 
sentially governed by the magnetic field B± . The strong 
damping of Rabi oscillations is caused by fast HH spin 
relaxation at B± = 0.8 Tw B^ 2 (first peak in Fig. [TJ. 
A slight increase in B± drastically affects the dynamics 
of the HH spin (dotted line in Fig. Rabi oscillations 
become stable due to much longer spin relaxation times 
and, at B± = 0.865 T« Bj_, the effective dipole moment 
dso ~ which substantially decreases the Rabi frequency 

In conclusion, we have introduced an efficient method 
for spin detection and manipulation of a HH in a QD. 
Furthermore, this method could be applied to determine 
important parameters of HHs in QDs (such as the effec- 
tive (/-factor, mass, SOI constants, spin relaxation and 
decoherence times) , from the position and shape of res- 
onances in the rf power absorbed by the system or from 



the amplitude evolution and the frequency of the Rabi 
oscillations. 
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